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Abstract. Using a mapping of compact polymers on the Manhattan lattice to 
spanning trees, we calculate exactly the average number of bends at infinite 
temperature. We then find, in a high temperature approximation, the energy of 
the system as a function of bending rigidity and polymer elasticity. We identify the 
universal mechanism for the relaxation of compact polymers and then endow the model 
with physically motivated dynamics in the convenient framework of the trees. We find 
aging and domain coarsening after quenches in temperature. We explain the slow 
dynamics in terms of the geometrical interconnections between the energy and the 
dynamics. 
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Introduction: Polymers confined to spaces much smaller than their radius of 
gyration are common in biological systems such as DNA packed in viral capsids p] or in 
the nucleus of a cell [2]. Both the geometrical structure and the dynamical behavior of 
strongly confined polymers are of considerable interest because of their role in biological 
functions. One of the first theoretical models for such phases is Flory's compact polymer 
model (CPM) [3] which describes a polymer by a Hamiltonian walk on a lattice. This 
two-dimensional model is relevant for polymers on surfaces [4] such as DNA adsorbed 
on a lipid bilayer. The thermodynamics of Flory compact polymers on the square lattice 
has been studied using field theoretical methods [5j [6] . Little is known, however, about 
the geometrical structure of compact polymers. Moreover, despite the fact that the 
Flory's CPM was introduced more than fifty years ago, the dynamics of the model has 
not been studied so far. 

In this letter, we use a mapping of the compact polymer on the Manhattan lattice 
to spanning trees [7J [8] to address the problems above. This mapping allows us to obtain 
an exact analytical expression for the average number of bends and then calculate the 
energy of the system in a high temperature approximation. More importantly, we 
uncover a deep universal connection between the dynamics of compact polymers and 
their geometric structure. We identify the extension and retraction of geometric features 
called fingers as the main mechanism for the relaxation of compact polymers. We then 
use the mapping to the spanning trees to endow the Flory's model with physically 
motivated kinetic rules. We find aging and slow relaxation in the quench dynamics of 
the system and argue that the slow relaxation stems from the geometrical correlations 
between the energy and the number of fingers. The relaxation phenomenology presented 
in this paper gives a universal description of polymer dynamics under strong geometric 
confinement. 

Model and Thermodynamics: The compact polymer is modeled as a Hamiltonian 
cycle, i.e. a self-avoiding loop which visits every site, on the Manhattan lattice (a 
square lattice with successive parallel lines oriented in alternating directions) with lattice 
spacing a/2. As seen in Fig. (TJ there is a one-to-one correspondence between polymer 
configurations and spanning trees on a square lattice with lattice spacing a. A spanning 
tree is a connected graph that visits every site of this new lattice without forming loops. 
Hereafter we work with the lattice underlying the trees and assume it has N = m x m 
sites. In a tree, a vertex can have coordination 1, 2, 3 and 4; sites with coordination 2 
can have the two bonds either "aligned" (coordination 2a) or "bent" (coordination 26). 
Denoting the number of sites with coordination a by N a , we have 

A x + A 2 + A 3 + A 4 = A 

Ax + 2A 2 + 3A 3 + 4A 4 = 2A (1) 

We assume periodic boundary conditions so there are 2A bonds on the lattice with 
(m — l)(m + 1) ~ A bonds belonging to the tree. For a given polymer configuration, 
we consider the bending energy and the elastic energy due to length fluctuations. As 
seen in Fig. El vertices of the tree with coordination 1, 2a, 26, 3 and 4 correspond to 



Letter to the Editor 



3 



t i 























i 


i 
> 














































1 

1 
























































-- 


i 






















ft 




1 










J L- 


— 1— - 









Figure 1. (Color online) A compact loop polymer (dashed red line) and the 
corresponding spanning tree (bold blue line). 



T 




Figure 2. (Color online) Configuration of the tree (bold blue line) and the polymer 
(dashed red line) at sites (black stars) with coordinations 1, 2a, 26, 3 and 4. 




Figure 3. Length fluctuation at a site with coordination 1. 



the polymer having 2, 0, 2, 2 and 4 bends respectively. Hence using Eq. ([I]), the total 
number of bends, is given by 

N B = 2(N + iV 4 - N 2a ) (2) 

If the characteristic energy of one bend is e&, the bending energy is given by Eb = SbNs- 
If we consider a polymer whose unstretched length is twice the length of the tree, and 
assume that the bends are rounded with a radius r, we can show using elementary 
geometry that the sites with coordination 1, 2a, 26, 3 and 4 contribute to the length 
fluctuation by a — 4r + nr, 0, — 4r + irr, —a — 4r + irr and —2a — 8r + 2-kt respectively. 
For example, the length fluctuation for a site with coordination 1 is equal to the length 
of the curve shown in Fig. [3] Using the above fluctuations and Eq. ([[]), it can be shown 
that the total length fluctuation depends only on N B and is given by |(7r — A)tNb- 
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By adding the corresponding elastic energy to the bending energy, we get for the total 
energy of the system 

E(N B ) = e b N B + ^(tt - m~) 2 N B 2 (3) 

where the spring constant of the polymer is ■Jfr-. To find the average number of bends 
(Nb), we need p 4 and p 2a where p a is the probability of a site belonging to the tree 
having coordination a . 

The probabilities pi, p 2 , P3 and P4 were calculated in reference [9] using Kirchhoff's 
matrix-tree theorem [JO]. Notice that because of constraints in Eq. ([TJ, only two of 
these four probabilities are independent. 

8 16 _ 8 20 16 

Pl=~2 3> P4=-l+ 2+^ ( 4 ) 

We also need p 2a , which we calculate below. Kirchhoff's theorem states that the number 
of spanning trees on a graph is given by any cofactor (they are all equal) of its Laplacian 
matrix. The Laplacian matrix of a graph with N vertices is an N by iV matrix whose 
diagonal elements are the coordination number of vertex % and the off-diagonal 
elements A^- are —1 times the number of links between vertices % and j. Representing 
empty, occupied and unconstrained bonds by , — and respectively, we can write 

p (.„{.„ ) -p{ A ) — p ( 1 )+p( [ ) (5) 

where p(- ■ ■) is the probability of the four bonds emanating from a site having the shown 
configuration. The first and second terms in the right hand side of Eq. [5] are each equal 
to pi/4 and the last term to p 2a /2. The left hand side of Eq. [5j p ( i ) = p c is equal to 
ratio of the number of spanning trees on a square lattice with two adjacent horizontal 
bonds removed and the total number of spanning trees on the original square lattice. 
Denoting the associated Laplacian matrices by A + B (the lattice minus two bonds) and 
A (the lattice), we have according to the matrix tree theorem p c = ^^j^p- Using the 
identity det(XY) = detXdetF, the ratio of cofactors can be written as det(J + GB) 
where G{r — r') = A _1 (r, r') — A~ 1 (r',r / ) is the regularized inverse Laplacian. Noticing 
that the Laplacian of the square lattice is diagonal in Fourier space we have [TT] 

G(±ax) = G(±ay) = -\ G(±ax ± ay) = -± 

G(±2ax) = G(±2ay) = -1 + | 

Since the matrix B describes a change in the connectivity of 3 sites, it only has a 
3x3 nonzero block. Then p c = det(J + GB) can be calculated by evaluating the 
determinant of a 3 x 3 matrix. We can then obtain p 2a using p\ from Eq. H] and the 
result is p 2a = 2p c — Pi — ^ — ^l + ^l, i n the limit of an infinite lattice. Thus in 
the thermodynamic limit, the average number of bends in a closed compact polymer 
(Nb) = 2N(1 + Pi — p 2a ) can be finally calculated and is given by 

{Nb)=n£-±). (7) 



(6) 
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The above probabilities have been computed assuming all configurations are equally 
likely (infinite temperature). For finite high temperatures the energy of the system can 
be obtained as follows: for a = 1,2a, 26, 3, 4, let us define f a (x) as the probability 
distribution function that a tree configuration has N a = x. Then f a (x) has a 
peak at Np a and we have numerically validated a Gaussian approximation for the 
fluctuations around the mean value of Np a with standard deviations scaling as a ~ y/~N 
by generating random spanning trees using Broders algorithm [12]. Similarly the 
distribution function of N B is approximated by a Gaussian peaked at (N B ) given by 

1 (iV B -<iVfl)) 2 

f(N B ) = ~^—e ™* (8) 
/2irNo-$ 



where <7q = 0.645 ± 0.003 has been obtained numerically. Within the Gaussian 
approximation (jSJ), the partition function Z oc ff^ e lN f(Ns)e~ l3E( ' NBS) dNB where the 
number of spanning trees on a square lattice with N sites asymptotically scales as 
e lN JT3] with 7 = — p + p — ^2 +■■■)■ From the partition function, the average energy 
of the system in the thermodynamic limit and at inverse temperature /3 is obtained as 

cvm AT (m - (3e b a^)(m e e + e b + (3e e e b a^) 

E{P) = N (l + 2^)> (9) 

where mo = ^jp- and e e = ff (tt— 4) 2 (^) 2 is a characteristic elastic energy. In the absence 
of length fluctuations, the energy reduces to E{j5) = NsbirriQ — fishO-^). Notice however 
that because of the high-temperature nature of the approximation, this equation does 
not hold for arbitrarily large flei, (the total energy cannot be negative) and it does not 
capture the melting transition the system [6]. 



The system has a melting transition [6] since the polymer is disordered at infinite 
temperature and the T = configurations are ordered as seen in Fig. HJ To find 
the transition temperature, we define the phenomenological order parameter Di = 
(N v — Nh)/(N V + Nh), with A^, and the number of vertical and horizontal bonds 
in a macroscopic (but not too large) region of size I. Calculating the Binder cumulant 
B e = 1 - pi] for I = 20, 28, 36, 44 to identify the transition we obtain C ~ 0.8. 

We calculated the cumulant using an ensemble average with 256 realizations obtained by 
Monte-Carlo (MC) simulations. Furthermore, the specific heat obtained numerically has 
a peak close to the transition whose finite size scaling indicates no thermodynamic limit 
divergence. This is suggestive of a Kosterlitz-Thouless transition [15]. The transition 
is analogous to that of [6] where the compact polymer lives on the square rather than 
the Manhattan lattice. 

Quench Dynamics and Aging: In this section we study the dynamics of a compact 
polymer. Notice that for a single strongly-confined polymer, the center of mass does 
not diffuse and the system dynamics is characterized only by reshaping. In order to 
study the dynamics of the system, we need to endow the Flory's model with some 
kinetic rules. These should represent the actual motion of the polymer rather than just 
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Figure 4. Some zero temperature configurations of the spanning tree which minimize 
the bending energy. 

transform the configurations as needed for thermodynamic MC studies. In reference [16J, 
we studied the dynamics of a confined single polymer using a microscopic model. We 
found that the reshaping of strongly confined polymers takes place primarily through 
fingering events in which a finger retracts and another finger extends to fill the void 
(Fig. [5]). In the language of the trees, a fingering event is a local bond flip at a site 




Figure 5. A fingering event. 



with coordination one: a bond is added to one of the three empty edges while the 
original bond emanating from the coordination-one site is removed from the tree. Since 
a fingering event involves the collective motion of many monomers, this dynamics is 
greatly accelerated. For a given move the change in the system energy AE depends 
only on the change in the coordination of the sites involved in the bond flip. To satisfy 
detailed balance we use Glauber-like dynamics where the rate r of a flip is proportional 
to e _/3A£; . As explained in reference |16j . although the shape of the polymer is modified 
locally, the actual polymer undergoes global motion through reptation [17] . Therefore, 
in addition to the energetics, a kinetic factor controls the rate at which the moves take 
place. To find the kinetic factor in the language of the trees, consider one fingering event. 
In a closed polymer there are two paths for monomers to flow from the retracting finger 
to the extending one; both paths contribute to the motion but in the thermodynamic 
limit, reptation essentially takes place along the shorter path of length i s . In order to 
create a new finger, the excess length coming from the retracting finger diffuses along 
the shorter path of length £ s with a diffusion constant proportional to j- [18J. Since 
the polymer reptates a fixed distance for each fingering event, we have r ex p Adding 
a bond to the tree (representing the extending finger) creates one and only one loop. 
This loop partitions the system into two regions of areas s = S and s = N — S, each 
containing one of the two polymer paths. For a path lying inside a region of area s, 
the length of the reptating polymer is then / = p — 2 + 2(s — p/2 + 1) = 2s where p 
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is the perimeter of the region and (s — p/2 + 1) is the number of lattice sites inside it. 
The rate r will then be proportional to m [ n ^g N s y The dynamics we have introduced 
above is "non-local" in the sense that for every move a non-local quantity S needs to 
be computed. A faster, "local" dynamics can be implemented by substituting - with 
the infinite-temperature average (-). The same thermodynamics and similar dynamical 
features as for the non-local dynamics are obtained. 

With open boundary conditions, we can show that, in the thermodynamic limit, 
the region with smaller area is the inside of the loop formed by adding a bond. At 
infinite temperature, this is a consequence of that fact the probability of forming loops 
of area s decays asymptotically as [5]. At low temperatures, the typical loop area 
5* scales as the linear size of the ordered domains which is much smaller than N — S. 

Using both the local and the non-local dynamics, we study the response of an ensemble 
of spanning trees to a sudden quench from infinite temperature to a non-vanishing 
inverse temperature 0. In the simulations we only consider the bending rigidity and 
neglect the elastic energy; temperatures are measured in natural units where £6 = 1. 
We measure how the average Nb oc E and iVi decay as a function of time. Since the 
dynamics takes place only through fingering events at the leaves of the tree (the mobile 
nodes), iVi is related to the system mobility. As seen in Fig. El for small values of /3 
the energy decays steadily to its equilibrium value whereas for larger a broad plateau 
emerges. N% has a similar behavior. All plateaus emerge at approximately the same 
value of ^? and for different temperatures. The results presented are for N = 28 2 . 
A parallel study for N = 18 2 indicates no system size dependence up to the time-scales 
we studied. For quenches to temperatures which allow the system to equilibrate within 
our simulation times, the system size independence of the relaxation time holds during 
the whole evolution. For quenches to lower temperatures, we observe that the time 
evolution is the same for the two different system sizes until the early stages of the 
secondary relaxation. The time scale necessary to reach complete ordering at such low 
temperatures is indeed expected to depend on the system size but such time scales are 
not accessible to our simulations. 

The dynamical features we observe can be explained by noting that none of the 
moves which decrease the energy increase the number of mobile nodes N\. Hence, as the 
system decreases its energy, the average distance between mobile nodes grows, freezing 
larger and larger parts of the system. Also, mobile nodes can be created or annihilated 
only if they are nearest neighbor to another mobile node and they can diffuse only by 
hopping over another mobile node. At zero temperature, only the moves that do not 
increase the energy are allowed and the plateau emerges once the moves that decrease 
the energy are exhausted. While in the plateau, the system is exploring its phase 
space with the moves that do not change the energy and a secondary decay begins only 
once the cooperative reshaping due to such moves has made further energy-decreasing 
moves accessible. At finite low temperatures the behavior is similar, although the finite 
probability of energy-increasing moves shortens the plateau. A relaxation mechanism 
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Figure 6. Emergence of a non-equilibrium plateau in the quench dynamics of the 
system. 

governed by fingers and the slow dynamics due to the geometric correlation between the 
number of fingers and the system energy are generic dynamical features which do not 
depend on the details of the model. Notably, with the non-local dynamics, iVg and Ni 
decay into a plateau of the same height as for the local dynamics. The decay, however, 
happens more slowly since the average loop area controlling the decay rate increases 
with time. 




Figure 7. Aging in the shape-shape correlation function c(t, t w ) for (3 = 6.5. 
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We now define a shape-shape correlation function for the system c(t,t w ) = 
-h X)f $i{t, tw), where Si(t,t w ) = 1 if the configuration of bonds emanating from site 
i are exactly the same at times t w and t and zero otherwise. As seen in Fig. [7J the 
system decorrelates more slowly for longer waiting times t w . The aging in the system 
is characterized by domain coarsening. As seen in Fig. El the ordered equilibrium state 
below the critical temperature consists primarily of a domain of horizontal or vertical 
bonds. When the system is driven out of equilibrium several ordered domains persist for 
long times (plateau) and undergo a slow coarsening during the secondary decay (aging). 




Figure 8. (Color online: thick black vertical bonds, thin green horizontal bonds) 
Left: a disordered equilibrium state after a quench from /3j = to /3f = 0.5 < j3 c . 
Center: an ordered equilibrium state for /3f — 1.0 > (3 C - Right: coexistence of domains 
in a non-equilibrium plateau for a quench to f3f = 3.0. 

Summary and Outlook: We used a mapping to spanning trees to calculate quantities 
related to the geometrical structure of Flory's compact polymers. We found that the 
dynamics of compact polymers is closely related to their geometric structure due to 
the role of fingers in reshaping. Formulating the dynamics in terms of the collective 
polymer reshaping, we endowed the model with physically motivated kinetics in the 
convenient framework of the trees. Our accelerated dynamics gives access to the long- 
time out-of-equilibrium phenomenology of compact polymers. We found that the system 
exhibits slow dynamics and aging characterized by the coarsening of ordered domains 
when quenched to low temperatures. We presented a generic phenomenological picture 
for the dynamics of strongly confined polymers by noting that the key to the slow 
relaxation is the geometrical correlations between the number of fingers and the system 
energy. Even in the continuum, where fingers do not have a rigorous definition, fingering 
events are a generic mechanism driving the relaxation dynamics of polymers [TSj and 
qualitative dynamical features such as domain coarsening and aging found in our paper 
are not expected to depend on the details of the lattice model. Quantitative results like 
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the number of bends however are specific to the lattice studied in the paper. Our model 
could be extended to study a fully-packed melt of loop polymers instead of single chain 
by using a mapping to a spanning forest and appropriate kinetic rules. Similar studies 
in the context of non-compact polymers have been done on a related problem of chains 
diffusing in an array of obstacles [191 EH El] . 

The simulations were carried out on Boston University supercomputer facilities. We 
are grateful to G. Biroli, C. Castelnovo, J. Jacobsen, J. Kondev, A. Meller, S. Redner, 
J. Schmit, and A. Sportiello for helpful discussions and correspondence. This work was 
supported in part by the NSF Grant DMR-0403997. 
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